At the very origin of chaos in three-dimensional dynamical systems 
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The mechanism responsible for the emergence of chaotic behavior has been identified analytically 
within a class of three-dimensional dynamical systems which generalize the well-known E.N. Lorenz 
1963 system. The dynamics in the phase space has been reformulated in terms of a first-exit-timc 
problem. Chaos emerges due to discontinuous solutions of a transcendental problem ruling the time 
for a particle to cross a potential wall. Numerical results point toward the genericity of the new 
found mechanism. 



PACS number(s) : 47.52. +j, 05.45.Ac 

Many chaotic phenomena in physical sciences and 
engineering can be satisfactorily described by low- 
dimensional autonomous dynamical systems. Among 
them the most celebrated example is the well-known 
E.N. Lorenz model |]|| 



a(y - x) 

—y + x + (r — 1)(1 — z)x 
~b(z - xy) , 



(1) 



here rewritten using rescaled variables x, y, z in the case 
r > 1. Besides it attracted the attention of the scientific 
community mainly on account of its ability to illustrate 
how a simple model can produce very rich and varied 
forms of dynamics, depending on the control parameters 
in the equations. In particular, chaotic behavior in its 
configuration space around three fixed unstable points 
(0, 0, 0) and (±1, ±1, 1) may take place for r larger than 
a critical value function of a and b [|],|| . 
Starting from the original work of Lorenz [jl| , a huge liter- 
ature has grown on the subject (see, e.g., Refs. for 
recent literature), almost independently from its some- 
what doubtful success of satisfactorily describing the 
fluid behavior related to the Rayleigh-Benard problem 
jjj. Nevertheless, for the class of Lorenz-like dynamical 
systems no determined attempt has been made yet in 
order to investigate analitycally the very origin of their 
chaotic behavior. This is the main aim of the present 
work. Specifically, we shall reformulate the Lorenz model 
dynamics in the phase space in terms of more familiar 
ideas such as mechanical properties of particles moving 
in one-dimensional potential fields subjected to viscous 
forcing. Such reformulation will be the starting point 
to generalize the original Lorenz system to a whole class 
of Lorenz-like systems. Among them we shall select a 
particular piecewise linearized model through which an- 
alytical results on the dynamics can be obtained. As we 
shall see, the dynamics generated by this model can be 
mapped into a first-exit-time problem. This will allow 
us to identify analytically the very origin of the system 
chaotic behavior. Finally, the robustness of the identi- 
fied mechanism will be tested numerically for other fully 
nonlinear generalized Lorenz-like systems. 



Let us start our analysis by noticing that, far from 
the initial transient, the system ([!]) is equivalent to the 
(integral-) differential equation 



' dx ■ 



(2) 



Eq. (0) can be interpreted as a customary classical dy- 
namical equation of motion for a (unit mass) particle 
subjected to a viscous force — (<r + l)x in the potential 
field U(x, t). Here, 



U(x,t) 



2a 



U(x) 



1 - 
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U t {x) 



(3) 



is a potential field resulting from a weighted average (in 
the chaotic regime b < 2a) of a constant quartic poten- 
tial U(x) = a(r — l)(x 2 — l) 2 /4 and a time depending 
quadratic potential U t (x) = a(r ~ l)[x 2 — l]b(x 2 — l)/2. 
We have defined 



[f] b = b dre-^f(t-r) 



(4) 



where f(t) is a generic function for which the integral 
exists. The latter integral clearly represents the memory 
of f(t) (at the time t), i.e. the present outcome of its ex- 
ponentially weighted past evolution. Thus, the potential 
Ut(x) depends on time through the exponential memory 
of the past motion. 

One can check that, given x from (|^), the variables y 
and z are obtained from the relations y = x/a + x and 
z = b/(2a)x 2 + (1 - b/(2a))[x 2 ]b, respectively. In this 
formulation it is for instance evident that x and y are 
synchronazing coordinates while this is not for z |^ . 

The above description immediately leads to the follow- 
ing considerations highlighting the role of the memory 
in the route to chaos. For Ut = 0, the particle motion 
stops after some time in one of the two minima of U{x). 
This because of the viscous term (a + l)x. On the con- 
trary, non trivial behaviors, including chaotic ones, may 
take place due to the statistical balance between dissipa- 
tion and energy exchanges produced by Ut- The bistable 
character of U plays a crucial role for the emergence of 
chaos. Indeed, initially very close trajectories starting in 



1 



the same cell, may undergo to completely different evo- 
lutions if, at a certain time, one has sufficiently energy 
to cross the peak in x = 0, while this is not for the other. 
This is the first clue that chaos arises from the unpre- 
dictability of the instants when particles pass through 
the barrier in x = 0. We shall give in the sequel the 
analytical proof of this heuristic argument together with 
the reason at the very origin of such unpredictability. To 
do that, let us give a further reformulation of the system 
( 1 ) . This will make it possible to generalize and simplify 
(1) with the final goal to deal with chaos analytically. 
By suitably scaling the time: t i— > r = [(r — f)/2] 1 / 2 t, 
Eq. (0) can be recast in the form M 



obtaining: 



x + rji + (x 2 — 1) x = —a [x 2 — f] 



(5) 



where rj, a and /3 are related to the original parameters 
<7, b and r. 

A more expressive form of Eq. (|^) is 



x + rjx + (q(x) + a [q{x)]^j $'(.t) = 



(6) 



where $>(x) = 1/2 a; 2 and q(x) 



•1 = 2$(x)-l. This 



equation describes the motion of a (unit mass) particle 
subjected to a viscous force — rjx and interacting with a 
fixed potential field <fr(x) through a 'dynamically varying 
charge' qt(x) — q{x) + a [q{x)] fj . It is constituted by a 
fixed 'core' charge, a 'locally acquired' charge, related to 
the local potential, and an exponentially vanishing 'mem- 
ory' charge, continuously depending on the previous in- 
stantaneous charge history. This scheme can be used 
to mimic the instantaneous effective charge of a particle 
moving in a background (structured) particle bath. In- 
deed, the coupling of [q(x)]^ with the background poten- 
tial <$>(x) yields an endogenous forcing term which allows 
a self-sustained unceasing motion, even in the presence 
of friction, and causes a corresponding unceasing inner 
transfer of the bath charge . 

At this point the Lorenz system ([!]) can be easily gen- 
eralized in the form 



x = a(y — x) 

y = —y + x+(r — l)(l — z)$'(x) 

z = -b [z - \q'{x){y - x) - q(x) - ll 



(7) 



from which Eq. (^|) follows after some algebraic manipu- 
lations. 

As already observed, the chaotic behavior of Lorenz 's 
system essentially depends on the unpredictability of the 
instants when x change its sign: as long as it keeps con- 
stant sign the system evolution is certainly nonlinear, and 
nevertheless not chaotic at all. This fact suggests a slight 
modification of the original form, in order to single out 
analytically the origin of chaos without to be faced with 
the difficulties arising from nonlinear problems. We thus 
set in Eqs. @ and (0) = \x\ and q(x) = $(ac) - I 



+ i]± + {\x\ - l + a [\x\ - 1]J sgn(a;) = (8) 



and the corresponding piecewise linearized system 
x = a(y — x) 

y = —y + x + (r — 1)(1 — z) sgn(a;) 
z = -bz + b sgn(x)^ 1 



(9) 



where sgn(x) = 

Our assumptions for $ and q correspond in Eq. (|^) to 
U(x) = air - \){\x\ - I) 2 /2 and U t (x) = a{r - l)[|z| - 
l]/3(|x| — I). The fundamental bistable character of U is 
thus maintained. 

In order to solve Eq. (||) , we exploit the fact that it is left 
invariant under the transformations x i— > —x, x i— > — x 
and x i— > —x. We can thus focus on one of the two 
regions x < and x > 0. Let us consider, e.g., x > . 
In this case Eq. (JsJ) is equivalent to the third order linear 
equation 



ft 
dr 3 



d 2 d 



di 



+ (/? + '?)7T + (1+^)t+^1 + «X = (10) 



dT 



dT 



where £ = x — 1. 

Notice that when the particles cross the barrier in x = 0, 
the evolution described by Eq. ( ^o|) has to restart with 
new initial conditions corresponding to an elastic collision 
against a rigid wall posed in x = 0. 

It can be checked that, in the case of chaotic motion we 
are interested in, the general solution of Eq. (10) can 
always be written in the form 

£(t) = e A " T (d cos(Air) + C 2 sin(A l r)) + C^- X ° T 

(11) 

with Ao,A r ,Ai > and C\, Ci and C3 real coefficients 
determined from the initial conditions on £, £ and £. 
The instant n at which the first particle collision against 
the wall in x — occurs is thus defined by the equation 
£(ti) = —1. From (f[l]), t\ is thus the smallest positive 
solution of the transcendental equation 



Cie A '' ri cos(A l T 1 ) + C* 2 e A '' T1 sin(A. J r 1 ) + 



C 3 e 



-Aoti 



-1. 



(12) 



Geometrically, we can interpret T\ as the first intersection 
of g{r) = -C 3 e~ XoT - 1 with ft(r) = C x e A '' r cos(Air) + 
C2 e XrT sin(Air). The function g is a decreasing exponen- 
tial and h an oscillating function with growing amplitude. 
It is thus easily understood why even a little modification 
of initial conditions can produce a discontinuous varia- 
tion of t\ (see Fig. 1). As we shall see, the same reason 
applies also for the class of systems @) . 
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FIG. 1. Graphical interpretation of the discontinuous char- 
acter of Ti for small changes of the initial conditions. Full 
line represents g(r), dashed lines h(r) for the parameters 
a — 6.50, /3 = 0.19, r\ — 0.78 and for initial conditions (a): 
Co = 2.26, £o = -2.00; (b): £ = 2.30, £ = -2.00. Bullets 
denote the first intersection between g and h defining the first 
collision time against the wall (a): ri = 3.91; (b) ri = 8.76. 

The dependence of n on the initial conditions Co and Co 
(Co = —1) is implicitly contained in Eq. (p^). Coefficients 
Ci, Ci and C3 are indeed linearly related to the initial 
conditions. The way to describe the global behavior of 
Ti in terms of the pair (Co, Co) although simple results 
quite lengthy. We thus confine our attention on the cor- 
responding graphical shape presented in Fig. 2. 




FIG. 2. The behavior of the first collision time, ri defined 
by Eq. (^), vs the initial conditions £0 and £0 for a = 6.50, 
P = 0.19, n = 0.78. 

From the figure it appears that t\ shows sensitivity with 
respect to the initial conditions only in a limited subset 
of the half-plane Co > 0. It is natural to expect that in 
the chaotic regime the system is quickly attracted inside 
this region. It is now clear why chaotic behaviors can 
arise also for apparently simple nonlinearities as isolated 
(non eliminable) discontinuities (see, e.g., Refs. 10-14|] 
for other examples of piecewise linearized systems dis- 
playing varied form of chaotic behaviors). 
Notice that the same behavior showed in Fig. 2 for t\ 



holds for the n-th collision time function of the 

(n — l)-th initial conditions (Co" >4o )■ 
For (Co, Co) = (Ao,— Aq) one has t\ — » 00 (the peak in 
Fig. 2). In this case, it is possible to show that the sys- 
tem configuration, described by C = (C( r ),C( r ),C( r )), 
exactly lies on the stable manifold W s = {t w 3 \ t > — 1} 
with 11)3 = (1, — Ao, Xq) and its motion is an exponential 
decay on the fixed point C = 0. 

For (Co, Co) (Ao, — Aq) the evolution in the phase space 
C, C, C consists of both a 'rapid' decay towards C = 
along the stable manifold W s and a 'slow' amplified os- 
cillation on the two-dimensional unstable manifold, W", 
generated by «q = (1, A c , A^) and W2 = (1, A c , A^), with 
A c = A r + i\i and A c = A r — iXi. 

Focusing now on the n-th collision, we thus have from 
Eq. ( p^[ ) the behavior of r„ vs the pair (Cg _1 , Co ~^)- 
Taking the first and the second time derivative of 
Eq. ([l2]), we can relate Co™^ an d Co™'' to the initial con- 
ditions Co" 1 ' ) an d Co™ 1 (remember that Ci, C2 and 
C3 are linearly related to such initial conditions). The 
result is a two-dimensional Poincare map connecting 
(Co >Co 1 ' ) ) ^0 (Co>Co )• ^he nonlinear character 
of the map is entirely contained in the transcendental 
and discontinuos relation between r n and the initial con- 
ditions. 
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FIG. 3. The one-dimensional map r n = T n (£ n V> ); (a): for 
the piecewise linearized system (^); (b): for the Lorenz model 



Let us now assume that the attraction towards the un- 
stable manifold W is sufficiently fast (the goodness of 
this approximation is controlled by (3 + rf) for the sys- 
tem to be considered as belonging to W u at the collision 
time against the plane C = — 1- If this is the case, 'im- 
mediately' after each collision the system thus necessary 
lies on a straight line. This implies a linear relation- 
ship between velocity and acceleration immediately after 
the n-th collision. The functional dependence of r n on 



so 



(n-l) 



and 



is thus restricted on a straight line in 



the plane Co 



(n-l) 



?0 



(n-l) 



This amounts to say that r n 

l(n-l) 



is a function of one variable alone, say, Co ■ The map 
T n = 7"n(Co ) is shown in Fig. 3(a). 
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Furthermore, the Poincare map con- 

necting (£o" _1 \ lei ^) t° (£o™\ Co™^ ) reduces to a one- 
dimensional map, say, between i^ 1 ^ an d £o™^- Fig- 4(a) 
shows this map for a = 6.50, (3 = 0.19, 77 = 0.78. Its 
behavior is evidently chaotic: the discontinuities derive 
from the analogous ones in the function r„ = t„(^q™ ). 
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FIG. 4. The one-dimensional map between £ and £q ~~ ; 
(a): for the piecewise linearized system (^|); (b): for the 
Lorenz model (|l|). 

The final question to be addressed concerns the rela- 
tion of our results with the original Lorenz system (|l|) 
and more generally with the class of dynamical systems 
(f^). One immediately realizes that, under the condition 
$ (0) = 0, the linear relationship between velocity and 
acceleration in x — follows: 



+ rjx { n) = 0. 



(13) 



As for the piecewise linearized system (|^) , the generalized 
Lorenz systems (fy are thus described by one-dimensional 
maps. Focusing in particular on the Lorenz model (Q), 
the analogous of the maps reported in Figs. 3(a) and 4(a) 
are shown in Figs. 3(b) and 4(b). Similar behaviors have 
been obtained for other choices of $ and g in (Q) . These 
maps have been derived by numerical integration of (Q) 
by a standard Runge-Kutta scheme, whereas we remem- 
ber that all results relative to the piecewise linearized 
system have been obtained analytically. 
The resemblance of Figs. 3(a), 4(a) with Figs. 3(b), 4(b) 
points toward the robustness of the mechanism we have 
identified as cause of chaos in the linearized system (^) . 

In conclusion, the very origin of chaos for a whole class 
of three-dimensional autonomous dynamical systems has 
been identified. Chaos is entirely contained in a transcen- 
dental equation ruling a first-exit-time problem whose 
solutions appear discontinuous for small changes in the 
initial conditions. Results have been obtained analyti- 
cally for a piecewise linearized model belonging to a more 
general class of dynamical systems. We however showed 
numerically that the basic reason for the chaos to emerge 
applies also for the general case. 

We thank G. Cassinelli, A. Celani, M. La Camera, 
M. Vergassola and A. Vulpiani for illuminating discus- 
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